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Abstract — This  work  focuses  on  the  study  of  the  flow  of 
air  in  natural  convection  in  a solar  tower  of  small  size. 
The  behavior  of  air  in  the  tower,  considered  as  a solar 
dryer,  provides  information  on  the  amount  of  heat 
absorbed  by  the  air  upon  entry  into  the  collector.  A 
theoretical  approach  allows  us  to  theoretically  simulate 
the  flow  by  using  a mathematical  model  characterizing 
the  physical  parameters  of  the  system  during  a daily 
sunshine.  An  analysis  of  this  phenomenon  is  made  and 
results  are  obtained. 

Keywords — Flow  of  air,  natural  convection,  solar 
tower,  dryer. 

I.  INTRODUCTION 

The  concept  of  solar  tower  or  solar  power  chimney  was 
proposed  first  by  Cabanyes  [1].  Solar  towers  consist  of  a 
collector,  a chimney  and  a turbine.  They  allow  the 
production  of  electricity.  A Solar  tower  without  its  turbine 
can  have  a use  other  than  that  of  power  generation:  the 
drying  of  food  products,  thanks  to  the  drying  air  coming 
out  of  its  collector;  hence  its  name  of  solar  tower  dryer. 
The  drying  of  food  products  is  an  appropriate  means  to 
curb  losses  of  food  after  harvest  and  expand  the 
consumption  of  these  products  during  periods  of  non- 
production. This  is  an  operation  which  is  much  practiced 
in  West  Africa  by  farmers.  Several  studies  have  shown 
the  difficulty  of  getting  indisputable  elements  of 
validation  and  comparison  between  different  solar  dryer 
models  due  to  its  technological  diversity. 

In  our  work,  we  study  the  operation  of  a solar  tower 
dryer,  running  empty.  The  objective  is  the  use  of  the 
energy  gain  received  by  the  air  through  the  collector  for 
the  drying  of  food  products.  To  do  so,  it  is  essential  to 
carry  out  a study  of  fluid  flow  in  the  system  in  order  to 
know  the  amount  of  energy  that  the  air  must  absorb  with 
an  average  humidity  of  40%  in  the  cases  of  BURKINA 
FASO. 


II.  DESCRIPTION  OF  THE  SOLAR  TOWER 

DRYER 

The  solar  tower  dryer  is  a natural  thermal  power 
generator,  which  uses  solar  radiation  to  increase  the 
internal  energy  of  the  air  flowing  in  its  collector.  The 
useful  gain  of  the  solar  collector  is  converted  into  kinetic 
energy  flow  circulating  in  the  chimney  where  trays  or 
racks  are  arranged. 

2.1.  Description  of  the  solar  tower 

As  part  of  our  work,  the  solar  tower  has  a function  other 
than  to  produce  electricity.  Used  as  solar  dryer,  it  is  small 
compared  to  that  producing  electricity  and  consists 
mainly  of  a collector  and  a chimney. 

> The  collector  has  a radius  of  1.2  meters.  This  is  the 
part  where  the  air  circulates.  It  is  heated  by 
absorption  of  solar  flux  and  subjected  to  phenomena 
of  convection.  The  collector  consists  of  a cover  of 
four  square  shaped  glass  windows  of  a square  meter 
and  an  absorber  in  alu-zinc  painted  matte  black, 
between  which  the  air  flows. 

r The  chimney  has  a diameter  of  0.4  meters  and  a 
height  of  3 meters;  it  is  the  drying  chamber  and  has 
racks  for  display  of  products  to  dry. 

The  chimney  has  two  roles:  that  of  drying-chamber 
and  that  of  drawing  of  drying  air.  The  greater  its 
height,  the  more  important  its  pull  effect  and  it 
promotes  the  renewal  of  air  in  the  [2-3]  system. 

2.2.  Functioning  of  the  solar  dryer 

The  products  to  be  dried  are  not  directly  exposed  to  solar 
radiation  and  that  makes  of  the  tower  an  indirect  solar 
dryer.  The  air  entering  the  collector  is  heated  by  natural 
convection  from  the  sensor  and  moves  through  the 
chimney  effect  towards  the  drying  chamber  where  the 
drying  air  takes  away  the  water  molecules  from  the 
products  to  be  dried  ( Fig.l). 
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Fig.l:  Diagram  of  the  experimental  solar  tower  dryer 


III.  MATHEMATICAL  MODELING  OF  THE 
AIR  FLOW  IN  THE  SOLAR  TOWER 
DRYER 

The  solar  dryer  is  said  to  be  empty  when  the  drying 
chamber  does  not  contain  any  products.  The  modeling 
consists  in  transcribing  physical  phenomena  that  can  be 
observed  and  equations  that  can  arise  at  different  levels  of 
the  system. 

3.1.  Simplifying  assumptions 

The  study  of  heat  transfer  in  the  system  can  be  simplified 
by  assuming: 

> A uni-dimensional  flow  depending  on  an  average 
radius  in  the  chimney  and  the  height  in  the  chimney. 

> Properties  of  the  air  such  as  dynamic  viscosity, 
thermal  conductivity,  density  and  heat  capacity  that 
are  considered  constant. 

> The  air  follows  the  Boussinesq  approximation  rules, 
which  results  in  the  expression: 

p = p0[l  - P(T  - T0)] 
with  (3  = — 

To 

> Transfer  by  conduction  of  the  sensor  (absorber)  is 
negligible  as  it  is  a thin  metal  and  its  conductivity 


coefficient  is  very  high;  the  temperatures  of  these 
two  faces  are  therefore  uniform. 

> The  heat  transfer  by  conduction  through  the  glass 
pane  is  negligible  because  it  has  a small  thickness; 
the  two  sides  have  substantially  the  same 
temperature  when  the  glass  is  exposed  to  sunlight. 

> The  radiative  exchanges  of  air  are  negligible. 


3.2.  Equation  governing  the  flow  of  air  in  the  system 
The  rectangular  collector  is  considered  circular  for  its 
modeling  in  the  system  of  polar  coordinates.  The 
equations  of  the  ah'  flow  are  a function  of  the  radial  axis  r. 
In  the  chimney,  the  equations  are  a function  of  the  axial 
axis  z. 

3.2.1.  Continuity  equation 

In  natural  convection,  air  velocity  is  ‘a  priori’ 
unknown.  So  we  assume  that  the  air  flow  is  constant  in 
the  system. 

V.S  = cte 

3.2.2.  Flow  equation 

dV  1 dP  T p 


— = J.g  - — g sin  a 

dt  p0  dr  p0 


(4) 
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Where 


J = 


f-v2 

2 D-g 


J is  the  pressure  drop  per  unit  length  and  f the  factor  of 
friction. 

3.2.3.  Heat  equation 


A heat  balance  of  the  collector  and  the  chimney  is 
required. 

> Heat  balance  for  the  collector 

Considering  the  above  assumptions  for  an 
elementary  control  volume  of  the  collector,  the 
following  Figure  (Fig. 2)  illustrates  the  mechanism  of 
heat  exchange. 


Tc 


Ar 

Fig.  2:  Illustration  diagram  of  the  different  heat  exchanges  in  the  collector 


N Thermal  balance  of  the  glass 

For  an  elementary  surface  and  space  A r directed 
towards  the  center  of  the  chimney,  we  get  the 
following  equation  (6): 

mvCpv5  = -SvhCvab(Tv  - Tab)  - Svhrvc(Tv  - 
Tc)  - SvUav(Tv  - Tab)  + SvhCvf(Tf  - Tv)  + 
Sphrpv(Tp  — Tv)  + SvavG 
N Thermal  balance  of  the  absorber 

The  solar  energy  received  by  the  absorber  is  the 
difference  between  that  transmitted  by  the  surface  of 


the  glass  and  that  transferred  to  the  fluid,  including 
losses. 

d'T 

mPCPP  if  = -Sphc_pf(Tp  - Tf)  - Sphr  pv(Tp  - 
Tv)  - SpUar(Tp  - Tab)  + SpapxvG 
N Thermal  balance  of  the  fluid  (air) 

The  heat  balance  of  the  fluid  inside  the  chimney  is 
given  by  the  following  equation: 

mf^Pf"it^  = ~Svhc_vf(Tf  - Tv)  + Sphc  pf(Tp  - Tf)  - 
mfSf(Tfe  - Tfs) 


> Heat  balance  of  the  drying  chamber  or  chimney 

The  various  modes  of  heat  exchange  for  a section  of  the  fireplace  are  show  in  (Fig.  3). 
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Air  outlet 


Fig.  3:  Illustration  diagram  of  the  different  thermal  exchanges  in  the  drying  chamber 


Thermal  balance  of  the  outer  wall  (pe) 

ttlnpCnnp  “ 5pehc  pe  ab  (Tpg  Tab) 


•pe'-ppe  gt 


Spehr_pec(Tpe  _ Tc)  + SpehCond(Tpi  ~ Tpe) 

Thermal  balance  of  the  inner  wall  (pi) 

piuc_pi  f(  H — lpiJ  ~ ,3pi11cond  (Tpi 


(9) 


dTp 

■pi'-ppi  at 


m„iC„r,i  ^ = S^h,  f(Tf  - Tpi)  - S„ih, 


Tpe) 

Thermal  balance  of  the  air  in  the  drying  chamber 

Heat  exchanges  related  to  the  air  are  described  by 
the  following  equation: 

9Tf 

dt 


mtCpf^  = -Spider, 


Tpi) 


nifSf(Tfe  - Tfs) 


3.3  Coefficients  of  thermal  exchanges 

3.3.1.  Coefficients  of  exchanges  by  convection 

The  coefficient  of  convective  heat  exchange  between 
the  environment  air  and  the  upper  wall  of  the  glass  is 
given  by  Me  Adams  equation  (1954)  [4]: 
hc_vab  = 5.67  + 3.86Vvent 

For  that  between  the  fluid  and  the  bottom  wall  of  the 
glass,  we  have  the  expression: 

hc  vf  = -r- 

For  that  between  the  absorber  and  the  fluid,  we  have 
the  same  expression  as  the  previous  one: 

hc-pf  - TT 

3.3.2.  Coefficients  of  exchange  by  radiation 
Coefficient  of  radiative  exchange  between  the  glass 
and  the  sky: 


hr_vc  - ct£v(Tv  + TC)(T2  + Tc2) 

With  Swinbank  relationship  [5] 

Tc  = 0,0552  Taj,5 

> Coefficient  of  radiative  exchange  between  the  glass 
and  the  absorber: 

u _ d(Tv+Tp)(T2+T2) 
nr_pv  — JL+JL-! 

£v  £p 

> Coefficient  of  radiative  exchange  between  the 
chimney  and  the  sky: 

hr_pec  = CT£pe(Tpe  + Tc)(Tpe  + Tc2)  (18) 

3.3.3.  Coefficients  of  loss 

> Loss  ratio  before  the  collector 

Uav  — hc  vab  T hr  vc 

> Coefficient  of  rear  losses  of  the  collector 

1 

Uav  — ep  t eisolant  , eplaque  , 1 (^0) 

^-p  Nsolant  ^-plaque  ^c_vab 

3.4.  Correlations  for  calculating  the  Nusselt  Number 
Nusselt  number  depends  on  the  Prandtl  number  and 
Grashof  number  or  Rayleigh  number  (Ra)  in  natural 
convection. 

The  absorber  and  the  glass  being  horizontal  plates  heated 
from  above,  the  Nusselt  number  of  equations  (11)  and 
(12)  is  determined  through  the  following  correlations  [6]. 
Nu  = 0.54.  Ra'25  (21) 

With  2.104  < Ra  < 8.106 

Nu  = 0.15Ra'33  (23) 

With  8.106  < Ra  < 8.1011  (24) 
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3.5.  Choice  of  the  method  of  resolution 
To  solve  the  equations  of  the  system,  we  made  the  choice 
of  the  method  of  numerical  solution  with  explicit  finite 
differences.  Using  the  method  by  slice,  the  system  is 
divided  into  several  dummy  slices  in  the  direction  of  the 
air  flow. 

IV.  RESULTS  AND  DISCUSSION 
4.1.  Materials  and  methods 


The  experiment  took  place  on  June  30,  2016  from  9 to  17 
o’clock.  We  used  a temperature  logger  named  midi-logger 
with  ten  lanes  equipped  with  Type  K thermocouples  for 
measuring  the  temperatures  of  the  glass,  absorber,  drying 
air  in  the  collector  and  the  chimney.  A radiometer  and  an 
anemometer  allowed  us  to  measure  respectively  the  solar 
radiation  of  the  day  and  the  air  speed  at  the  exit  of  the 
chimney.  The  figure  (Fig. 4)  below  shows  a photo  of  the 
experimental  device. 


Air  outlet 


Chimney 
or  drying 
chamber 


Air  inlet 


Fig. 4:  Photo  of  the  experimental  device 


4.2.  Comparison  of  theoretical  and  experimental  global 
solar  radiation 

The  following  figures  (Fig. 5)  and  (Fig. 6)  respectively 
show  the  variation  of  the  theoretical  and  experimental 
global  solar  radiation  depending  on  the  time  of  day  on  30 
June  2016  in  Ouagadougou. 


Fig. 5:  Theoretical  variation  in  solar  radiation 
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Global  radiation 


Fig.  6:  Experimental  Variation  of  solar  radiation 
The  manipulation  began  around  09  o’clock  and  we 
already  saw  a flux  of  650W/m2  slightly  lower  than  that 
theoretically  expected  to  be  700W/m2  at  the  same  time. 
Then  this  flow  increases  to  reach  its  peak  of  950W/m2  at 
noon.  We  observe  a sinusoidal  growth  in  theory  but  the 
peak  at  12o’clock  is  967W/m2.  Some  cloud  passages  were 
observed  in  the  said  day.  The  average  systematic  error 
(RMSE)  between  the  theoretical  and  experimental  values 
is  about  56  W/m2. 

4.3.  Experimental  and  simulation  results  of  the  collector 
4.3.1.  Results  presentations 

Figures  (Fig. 7)  and  (Fig. 8)  show  the  variation  of  the 
average  temperatures  of  the  glass,  air  and  absorber  as  a 
function  of  time. 


Fig.  7:  Theoretical  temperature  of  the  air,  the  glass  and 
the  absorber  in  the  collector 


Fig. 8:  Experimental  temperatures  of  air,  glass  and 
absorber  in  the  collector 

The  various  average  temperatures  of  the  absorber  are 
superior  to  those  of  the  air  and  glass  respectively.  We  find 
that  the  temperatures  rise  gradually  as  solar  radiation 
increases  to  maximum  values  between  12o’clock  and 
13o’clock.  The  theoretical  maximum  temperature  of  the 
absorber  is  81°  C against  88°  C in  practice.  That  of  the 
glass  is  49°  C against  55  ° C and  that  of  the  air  63  ° C 
against  62°  C. 

4.3.2)  Comparison  between  the  various  temperatures  of 
the  air  at  the  outlet  of  the  collector 
To  validate  our  results,  a comparative  study  was  carried 
out  with  experimental  and  theoretical  results.  The  figure 
(Fig. 9)  below  shows  us  that  comparison. 


Experimental  temperatures  of  the  drying  air 
Theoretical  Temperatures  of  the  drying  air 


Fig. 9:  Comparison  of  experimental  and  theoretical 
temperature  of  the  drying  air 
From  9 o’clock  to  13.50,  we  see  that  the  predicted 
temperatures  are  higher  than  the  experimental 
temperatures  while  from  17o’clock  to  14o’clock, 
experimental  temperatures  are  slightly  higher  than  those 
of  the  simulation.  They  increase  or  decrease  according  to 
the  sunlight  of  that  day  and  reach  a maximum  value  of 
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about  62  ° C for  the  experimental  temperature,  and  about 
63  ° C for  the  simulated  one  with  RMSE  = 5.85  ° C.  This 
same  observation  was  made  by  al  Hakim  and  sowed  in 
2013,  when  he  simulated  a solar  tower  with  and  without 
storage  system.  [7] 

The  temperature  of  drying  air  obtained  at  the  outlet  of  the 
chimney  enables  the  drying  of  certain  food  products 
whose  drying  temperatures  do  not  exceed  62  ° C. 

4.4.  Temperature  Fields  at  the  level  of  the  racks 
The  average  temperatures  of  the  drying  air  at  the  racks  are 
given  in  Figures  (Fig. 10)  and  (Fig.l  1).  They  respectively 
show  theoretical  and  experimental  temperatures  at  the 
four  racks  in  function  of  time. 


Fig. 10:  Theoretical  average  temperature  of  hot  air  in 
racks 


Time(h) 

Fig.l  1:  Experimental  mean  temperature  of  hot  air  at  the 
racks 

The  temperature  of  the  drying  air  at  the  racks  also 
depends  on  sunshine.  Rack  N°  1 is  traversed  by  a flow  of 
hot  air  of  higher  temperature  compared  to  the  others.  We 


ISSN:  2349-6495(P)  I 2456-1908(0) 
notice  a decrease  in  the  temperature  of  the  hot  air  at  the 
level  of  the  racks  (from  rack  N°  1 to  N°  4).  Respectively 
it  shows  theoretical  and  experimental  maximum 
temperature  of: 

> 58  ° C and  59  ° C for  rack  N°  1, 

> 55  ° C and  56  ° C for  rack  N°  2, 

> 51  ° C and  52  ° C for  rack  N°  3 

> 47  ° C and  48  ° C for  rack  N°  4. 

This  decrease  in  air  temperature  was  noted  by  O. 
MOCTAR  et  al  in  2015,  when  he  stated  that  the  air 
temperature  in  the  chimney  decreases  gradually  as  the  air 
is  ascending. 

This  range  of  air  temperature  at  the  location  of  the  racks 
will  allow  us  to  dry  products  such  as  okra,  medicinal  or 
flavor  leaves,  green  peas  etc.  [5] 

4.5.  Velocity  field  outlet  of  the  chimney 
In  the  design  of  solar  dryer  dimensions,  the  temperature 
of  the  drying  air  and  humidity  are  certainly  important,  but 
so  does  the  velocity  of  the  air;  hence  the  need  to  compare 
the  simulated  velocity  to  the  experimental  ones.  Figure 
(Fig.  12)  shows  us  the  comparison  of  theoretical  and 
experimental  speeds  out  of  the  chimney  as  a function  of 
time. 


Experimental  air  velocity  at  the  chimney  outlet 
Theoretical  air  velocity  at  the  chimney  outlet 


Fig.  12:  Comparison  of  experimental  and  theoretical 
speeds  at  the  exit  of  the  chimney. 

The  fluctuations  of  experimental  velocity  are  due  to  the 
observed  passages  of  cloud  because  these  velocity  also 
depend  on  sunlight.  With  a mean  bias  of  0.18  m / s,  the 
maximum  value  of  velocity  measured  at  the  outlet  of  the 
chimney  or  the  drying  chamber  is  0.81m  / s while  that  of 
the  theoretical  is  in  the  range  of  0.72m  Is. 

V.  CONCLUSION 

Knowledge  of  the  behavior  of  the  airflow  in  the  solar 
tower  dryer  is  very  important  for  the  drying  process.  The 
range  of  theoretical  and  experimental  temperatures  at 
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different  racks,  tells  us  about  the  type  of  products  to  dry. 
Experimental  velocity  (0.81m  / s)  and  theoretical  (0.72m  / 
s)  out  of  the  room  to  dry  are  low  because  our  system 
works  by  natural  convection.  They  have  an  impact  on  the 
drying  time  because  they  do  not  promote  quick  drying  but 
can  prevent  the  crusting  of  the  product.  These  results  will 
test  the  solar  tower  dryer  by  performing  the  drying  of 
some  agricultural  products  and  improve  our  IT  program 
to  reduce  average  systematic  errors. 

Nomenclature 
S:  Surface  (m2) 

Cp:  Specific  heat  at  constant  pressure  (kj.  kg-1K-1) 
g:  Acceleration  of  gravity  (m.  s_1) 

G:  Golar  radiation  (m.  s_1) 
m:  Mass  (kg) 

h:  Heat  transfer  coefficient  W.  m-2) 

' m:  Mass  flow  ( kg.  s-1) 

R:  radius  of  collector  m 
P:  pressure  (Pa) 

Nu:  Nusselt  Number 
T:  Temperature  (K) 

V:  velocity  ( m.  s-1) 
dr:  Variable  Radius 
dz:  Variable  height 
e:  Thickness  (m) 

U:  Loss  Factor  W.  m-2) 

J:  losses 
f:  friction  factor 
Greek  symbols 

k:  Thermal  conductivity  (W.  m-1K-1) 
x:  Transmissivity 
e:  Emissivity 

p:  Density  of  the  air  mass  (kg.  m~3) 


a:  Stefan-Boltzmann  constant  W.  m-2k-4) 

P:  Coefficient  of  thermal  expansion  (K-1) 

Indices 
ab:  Ambient 
p:  Absorber 
f:  fluid  (air) 
v:  Glass  window 
av:  Front  window 
RMSE:  average  bias 

REFERENCES 

[1]  I.  Cabanyes,  "Proyecto  de  motor  solar,  la  energia 
electrica ",  Revista  general  de  electricidad  y sus 
aplicaciones,  Vol  8,pp.  61-65,  1903. 

[2]  S.  Roozbeh,  A.  Majid,  H.  Behzad,  ” Modeling  and 
numerical  simulation  of  solar  chimney  power 
plants”,  Solar  Energy,  Vol.  85,  pp.  829-838  , May 
2011. 

[3]  M.  Ousmane,  B.  Dianda,  S.  Kam,  “Experimental 
study  in  natural  convection”,  Global  journal  of  pure 
and  applied  sciences,  Vol  21, pp.  155-169,  2015 

[4]  M.  Daguenet,  Les  sechoirs  solaires  : theorie  et 
pratique.  Unesco,  1985. 

[5]  K.  Swinbank.  “Thermal  performance  of  solar  air 
heater:  Mathematical  model  and  solution  procedure 
solar  energy",  vol  55,  1995. 

[6]  A.  Yunus  Cengel,  J.Afshin  Ghajar,  "Heat  and  Mass 
Transfer:  Fundamentals  & Applications."  McGraw- 
Hill,  2011,  p.  567. 

[7]  H.  Semai  and  A.  Bouhdjar , " Modelisation  d’une 
centrale  a cheminee  solaire  en  regime  turbulent  et 
avec  stockage  thermique  ",  16eme  journee 
internationale  de  thermique. 


www.ijaers.com 


Page  | 118 


